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«j | Two-way selections play an important role in human social life. In this paper, we analyze the minimum model 

of random two-way selections system and provide its exact results on the matching ratio. Defining two basic 
relative ratios: n for the ratio of selection quantities from both sides, and £ for the ratio of the total number of 
states to the smaller number of two sides, we find the matching rate P of two-way selections trends toward an 
inverse proportion to both rj and f , in agreement with the numerical simulations. We also surprisingly notice that, 
most of the empirical data of the matchmaking fairs, the typical real-world two-way selections, are completely 
included in the range predicted by our model. More interestingly, by the fitting of empirical data, we find that the 
average number of characters in real- world human sexual chosen is about 15, which is generally in agreement 
with daily experiences, implying our model would be useful for the estimation on the main factors of affecting 
the results of two-way selections. 
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^ ' I. INTRODUCTION 

Two-way selections between people are very common in daily life, such as the marriage selections between men and women, 
the selections between seekers and recruiters, and trades between buyers and sellers. In a sense, two-way selections can be 
<"*> |' regarded as the base in the build of many social relationships. Generally, the participants in a two-way selection process usually 
l— ~" ', consider many characters of the people in the other side. For instance, in the human marriage selections, her/his looks, person- 
ality, wealth, sense of humor, and so on, are usually in the considerations. These characters, which generally can not clearly be 
distinguished and are difficult to be estimated through traditional methods such as psychological test and social survey, deeply 
£f~) • affect the result of selection. Furthermore, two-way selections are usually affected by many factors, such as, the size of members 
on the two sides and the corresponding ratio. For example, the supply-demand ratio on market would deeply influence on the 
prices and business practices. Therefore, it would be significant to study the basic laws dominating the two-way selections, 
including the factors affecting the matching rate and the method for estimating the total number of characters affecting the se- 
lection. The basic laws may be useful for several social services such as the prediction of friendships and social relations|Q]-[9l], 
CO ■ seeking marital partners, and so on. 

In this paper, we propose a minimum model for two-way selections to investigate the factors influencing on the matching rate. 
Based on our model, we also find an efficient way to estimate the total number of characters affecting the two-way selections 
from empirical data. 



X 



n. THE MODEL AND RESULTS 

Generally speaking, two-way selections have at least two parts of members, and the condition of successful matching is 
usually that the two sides of pairing satisfy the characters for each other. The minimum model of random two-way selections 
can be described as follows: i). The systems has two types of agents, A and B, respectively with the total number k\ and Ic2. ii). 
Each agent (the z'-th agent, say), either in type A or type B, has a own character signed by Cm (z e A) or Csi (z e B). At the same 
time, the character which each agent want to select is signed by Sai or S Bi- Therefore, describing the state of each agent Az or Bi 
requires two characters which are independently chosen from n kinds of characters following the uniform distribution, iii). The 
condition of successful match of two agents (z € A and j e B) is Cm = S Bj and Sm = Csj- 
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For given ki, ki, and n, the total number of matching pairs E can be written out according to the rules above: 
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According to : 

1=0 v ' i=0 

where /I = x£, Bessell(0, 2 -0c) and Bessell(l, 2 y[x) are Bessel functions of the amount of virtual cases. From Eq.(Q]i we obtain 
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where /li = k\/n 2 , Xq. = A^/n 2 . 

Due to the symmetry of kl and k2 in Eq.(l) , it is enough to study the range k\ < ki- We discuss Eq.([T]i under three conditions 
of different relative magnitude of k\,ko and n 2 . When k\ < ki <sc n 2 , resulting in A[ — > 0, A2 — > 0, we only need to calculate the 
zeroth power term and the first power term of Eq.(|3]. We obtain 
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E =s — rfclfe. 



When < n 2 «; &2 or n 2 « < £2, according to 
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lim Bessell(0, 2jc) = lim Bessell(l, 2x) = 
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from Eq.© we obtain 



E » ki - k\ 
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Because Eq.© is obtained under the condition of £1 < « 2 <k £2 or n 2 « < &2, ^2 is very large, from Eq.(fj) we can obtain 
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We defined that 
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Therefore, 77 represents the ratio between £2 and k\ ; £ represents the ratio between n 2 and k t ; P represents the probability of the 
agents which successfully match. The Eq.© is transformed into 



P = 
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The Eq.Q can written as: 



The Eq.(0 can written as: 
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(1 < £ <sc 77 or £ <sc 1 < 77). 
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FIG. 1: (Color online) The comparison between analytical predictions and the simulation results. In the two sub-figures, the parameter 
ki = 100; rj is assigned values from 1 to 1000; £ is assigned values from 0.1 to 1000. (a) shows analytical predictions of Eq.lQJ; (b) shows the 
simulation results. 




FIG. 2: (Color online) The comparison between analytical predictions and the simulation results in the log-log plots. In the four sub-figures, 
the parameter fc, = 100; the squares indicate the simulation data; the solid lines indicate analytical predictions, (a), rj = 2; i; is assigned values 
from 100 to 1000; the solid line is obtained from Eq. J 10t . (b), £ = 2; r\ is assigned values from 10 to 1000; the solid line is obtained from 
Eq.QTJ. (c), rj = 1 ; tf is assigned values from 10 to 100; the solid line is obtained from Eq.dl2t. (d), f = 0.1; 77 is assigned values from 10 to 
1000; the solid line is obtained from Eq.{TQ. 



Fig l.(a) and Fig l.(b) show the comparison between analytical predictions of Eq.(Q} and the simulation results. Fig 2. (a) 
shows the comparison between analytical predictions of Eq.dTOb and the simulation results under condition 1 < 77 <K £ and 
displays a power-law distribution with the exponent -1 between P and Fig 2.(b) shows the comparison between analytical 
predictions of Eq. ([m > and the simulation results under condition 1 < £ <K 77 and displays a power-law distribution with the 
exponent -1 between P and 77; Fig 2.(d) shows the comparison between analytical predictions of Eq.dTTI) and the simulation 
results under condition £ <K 1 < 77 and displays the same power-law distribution to Fig 2.b. The above analytical predictions and 
simulation results are consistent with each other. It is that to say all analytical results are reliable. 

Consider a special case k\ — Ic2 — k, resulting in 77 = 1. On the one hand, from Eq.dTOt we obtain 

1 

P«-. (12) 
The relation between P and £ approximate to power law with the exponent -l,and this case is well shown in Fig 2.(c); from 
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FIG. 3: (Color online) The relationship between the experimental evidence and analytical predictions in the log-log plots. The red curve and 
the olive curve are obtained from Eq.l(9), and The parameters of red curve are r\ = 1, n = 15; The parameters of olive curve are rj = 356.67, 
n = 15. The round dots represent the empirical data in APPENDIX A. The £ mnv represent the maximum value 225 of 



Eq.dTTT) we obtain P m 1, it means that almost all of agents can match successfully under condition £ <sc 1 < 77. On the other 
hand, according to k\ -k^- k, resulting in Ai = A2, from Eq.© we obtain 



E~k- J— yfk. (13) 
\ n 

The second term of Eq.([T3l is the number of the agents which can not successfully match in type A or type B. The larger k is, 
The smaller ^frP~f^ Vk/k is. In reality, It's a result of fluctuation. There are n 2 kinds of the states that a agent maybe have in the 
model. In theory, the expected times of each state appearing is kin 1 . However, due to the fluctuations, almost all frequencies 
of every states appearing deviate around k/n 2 . As a result, some agents can not successfully match. The number of times that 
each state may appear obeys the binomial distribution. The fluctuation is close related to the standard deviation, according to the 
binomial theorem and standard deviation formula, we can obtain that the standard deviation equals to ^jk(\ - l/n 2 )/n 2 , which 
is directly proportional to Vk. Thus, the number of the agents which can not successfully match is also proportional to Vk. It 
explains the relationship between the second item of Eq.(fT3ll and Vic. From Eq.dOll. we know that the proportionality coefficient 
is ^n 2 In. 

Finally, we validate our model in the experimental evidence way. The marriage selections between men and women is a typical 
real-world two-way selections system. Here we collect eighty-two news of real-world matchmaking fairs from the Internet (see 
APPENDIX A). But the information of our data are incomplete. Only the numbers of participants and matching pairs are 
recorded. The numbers of men and women joined are omitted. So parameter 77 can not be determined exactly, but can only be 
restricted at a range. And n is a internal parameter need to be measured. We deal with the data and measure n below. 

We first consider the parameter 77. Because £1 < k%, the minimum value of 77 is n m j„ = 1 . For each set of data, there is a maximal 
value. Take the first set of data in TABLE I for example, there are 13 participants and 3 couples matched, thus k\ t i > 3, according 
to Eq.©, 77 < 77^1 = (13 - 3)/3. And we define two signs: r}maxjnm = minf^ T] max j = 2, r] max , max = max^ 77^^ = 365.67 by 
calculating eighty-two sets of data in TABLE I. We estimate n by minimizing 
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(PdataJ Panalysis,i) (14) 



in which Pdata comes from the experimental data and P analysis calculated from Eq|9] Firstly, we take 77 = 77„„„ = 1, and n 
ranges from 1 to 10 4 with a step size of 1, by calculating Eq.(fl~4b. we find that when 77 = 17, Eq.(fl~4b is minimal. Then 
we take 77 = r\ max , m in = 2, the same to way, we minimize Eq.(fT4b and find 77 = 13. In the Eq.©, n 2 and 77 are inversely 
proportional relationship , therefore, we define 7? mflA = 17 and 77„„„ = 13, and 77 belongs to [13, 17]. Finally we estimate n by 
n = (n max + n min )/2 = 15 and fix it. 
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Fig 3 shows the relationship between the experimental evidence and analytical predictions of our model. The red curve 
and olive curve are obtained from Eq.®. The parameters of red curve are rj — 1, n = 15; the parameters of olive curve are 
T] - 356.667, n - 15. From the Eq.©, when k\ is equal to the minimum 1, the £ take the maximum value 225. The error bar 
of ordinate P of round dots represents the range of empirical data P in TABLE I. Because k\ is unknown and ^ is undetermined, 
we can only obtain fci m ,„,; = E; < k\j < k\ max j = {k\ + ^)/2 and l5 2 /ki t „ tax j < £ < 15 2 /£i,mm,(- Therefore, the rang of abscissa 
£ of round dots is relatively wide and the middle points take = (15 2 /ki imax j + 15 2 /^i, mi >, ii )/2. 

We can see from Fig. 3 that when £ is relatively small, corresponding k\ relatively big, all empirical data are enclosed in 
between the the two curves; when £ is relatively big, corresponding k\ relatively small, some empirical data are enclosed in 
between the two curves, however others empirical data lie the red curve above and the trend of the empirical data is opposite to 
the analytical predictions. We have two reasons to explain them: one the one hand, organizers of some matchmaking fairs select 
less guests who meet their requirements to participate the matchmaking fairs from a lot of people who enter name and want to 
participate the matchmaking fairs. In this case, each guest is easier to find the right opposite sex who can match for each other; 
one the other hand, when the number of guests is small in a matchmaking fair, they know that it is very difficult to find the right 
opposite sex who can match for each other, so they are not to find the right opposite sex but tend to choose a basically mutually 
acceptable opposite sex. In this psychological effect, each guest is also easier to find the opposite sex to pair. The two reasons 
above cause that the fewer the guests are, the higher the probability p is. From the above analysis, the experimental data can be 
well explained by the model. 

III. CONCLUSIONS 

In summary, we provide the exact solution for random two-way selections system. The dynamics mechanism of a model 
is established and the analytical result is obtained. In several different conditions, we obtain the approximative and compact 
analytical results which agree with the simulation results. In the model, the parameter n is the total number of the characters 
from which each agent can chose and directly determine the probability of successful match. Therefore, n is a very important 
parameter for random two-way selections system. We propose a method to estimate the value of n by fitting empirical data. In 
this paper, the typical experimental data is collected by the Internet and the range of parameter n is estimated by the method we 
proposed. We find that the average of n is 15, which is in agreement with daily experiences, and that most of the experimental 
evidence data fall into the range predicted by our model. So the method would be useful to estimate the parameter n, and such 
simple model indeed catch the core of random two-way selections system in reality and would be valuable to the real life and 
several online social services. 
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Appendix A: Empirical Data of Matchmaking Fairs 

We collect 82 online news of real-world matchmaking fairs and record the number of guests and matching pairs. Due to some 
of news did not report exact numbers, we suppose each of these uncertain numbers has 10% possible ranges. The reports of 
these uncertain numbers generally have three types: nearly x, about x, and over x (here x is the number). For '"nearly x", we set 
the possible range is 0.95x + 0.05x, and 1 .00x ± 0.05* for "about x", and 1 .05 x ± 0.05x for "over x". All of these information is 
shown in Table 1 . 
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TABLE I: The data of matchmaking fairs 
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13 joined; 3 matched 
18 joined; 6 matched 

20 joined; 1 matched 

21 joined; 6 matched 

25 joined; 5 matched 

26 joined; 3 matched 
26 joined; 4 matched 
30 joined; 2 matched 
32 joined; 4 matched 
36 joined; 6 matched 
36 joined; 6 matched 
38 joined; 6 matched 
40 joined; 8 matched 

over 40 joined; 3 matched 
over 50 joined; 3 maiched 
aboul 60 joined; 5 maiched 
over 60 joined; 4 maiched 
over 60 joined; 5 maiched 
over 60 joined: aboul 10 maiched 
72 joined; 1 1 matched 
80 joined; 5 matched 
80joined; 10 maiched 
80 joined; 13 maiched 
80joined; 18 maiched 
nearly 100 joined; 5 maiched 
99 joined; 8 maiched 
IOOjoined; 5 maiched 
about IOOjoined; 7 maiched 
over IOOjoined; 16 maiched 

150 joined; 7 maiched 
nearly 200 joined: 8 maiched 
nearly 200 joined; 22 maiched 
about 200 joined; 4 maiched 

206 joined; 10 matched 
over 200 joined; 7 maiched 
over 200joined; 8 maiched 
over 200 joined: 38 maiched 

216joined; 19matched 
over 240 joined; 22 maiched 
about 258 joined; over 10 matched 
nearly 300 joined: 4 matched 
aboul 300 joined; 8 matched 
over 300 joined; over 10 matched 
over 300 joined: 32 maiched 
400 joined ; nearly 20 maiched 
over 500 joined; 3 maiched 
over 500 joined; 8 maiched 
over 500 joined; over 10 matched 
nearly 600 joined; nearly 40 matched 
over 600 joined; over 78 maiched 

nearly 800 joined; 58 matched 
over 800 joined; over 20 matched 
nearly 1000 joined; aboul 20 maiched 
nearly lOOOjoined; 58 matched 
nearly lOOOjoined; 64 maiched 
aboul lOOO joined; 12 matched 
about lOOOjoined; 15 maiched 
aboul 1000 joined: nearly 100 maiched 
over lOOOjoined; 3 maiched 
over 1000 joined; 4 maiched 
over 1500 joined; 48 maiched 
over ISOOjoined; over 100 maiched 

over 1600 joined; 31 matched 
over 2000 joined; nearly 100 maiched 
over 2000 joined; over 113 maiched 
nearly 3000 joined: nearly 100 maiched 
aboul 3000 joined; 186 
over 3000 joined; over 200 maiched 
over 4000 joined; over 500 maiched 
nearly 5000 joined; 108 maiched 
over 5000 joined; 218 maiched 
over 5000 joined; 231 matched 
over 5000 joined; 237 maiched 
over 6000 joined; over 270 maiched 
nearly lOOOOjoined; 28 maiched 
nearly 10000 joined; about 100 maiched 
nearly lOOOO joined; nearly 2000 maiched 
over lOOOOjoined; about 400 matched 
over lOOOO joined; over 1000 matched 
about 16000 joined; over 700 maiched 
over 16000 joined; over 600 maiched 
over 50000 joined; over 3000 matched 



42 ±2 
53+3 
60 ±3 
63 + 3 
63 + 3 
63 ±3 
72 



100 + 5 
110 ± 10 

150 
190 ± 10 
190 ± 10 
200 ± 10 

206 
210± 10 
210 ± 10 
210± 10 

216 
252 ± 12 
258 ± 13 
285 ± 15 
300 ± 15 
315+15 
315 + 15 

400 
525 ±25 
525 ± 25 
525 ± 25 
570 ± 30 
630 ± 30 
760 ± 40 
840 ± 40 
950 ± 50 
950 ± 50 
950 ± 50 
1000 ±50 
1000 ± 50 
1000 ±50 
1050 ± 50 
1050 ± 50 
1575 ±75 
1575 ± 75 
1680 ±80 
2100 ± 100 
2100 ± 100 
2850 ± 150 
3000 ± 150 
3150 ± 150 
4200 ± 200 
4750 ± 250 
5250 ± 250 
5250 ± 250 
5250 ± 250 
6300 ± 300 
9500 ± 500 
9500 ± 500 
9500 ± 500 
10500 ± 500 
10500 ±500 
16000 ± 800 
16800 ±800 
52500 ± 2500 



11 ± 1 

32 
19 ±1 



284 ± 14 

28 
100 ±5 
1900 ± 100 
400 ±20 
1050 ±50 
735 ± 35 
630 ± 30 
3150 ± 150 



0.4615 

0.6667 

0.1000 

0.5714 

0.4000 

0.2308 

0.3077 

0.1333 

0.2500 

0.3333 

0.3333 

0.3158 

0.4000 
0.1432 ±0.0068 
0.1 136 ±0.0064 
0.1671 ±0.0084 
0.1273 ±0.0061 
0.1591 ±0.0076 
0.3197 ±0.0470 

0.3056 

0.1250 

0.2500 

0.3250 

0.4500 
0.1056 ±0.0056 

0.1616 

0.1000 
0.1404 ±0.0070 
0.2933 ±0.0267 

0.0933 
0.0844 ± 0.0044 
0.2322 ± 0.0122 
0.0401 ±0.0020 

0.0971 
0.0668 ± 0.0032 
0.0764 ± 0.0036 
0.3627 ± 0.0173 

0.1759 
0.1750 ±0.0083 
0.0859 ±0.0121 
0.0282 ±0.0014 
0.0535 ± 0.0027 
0.0703 ± 0.0097 
0.2036 ± 0.0097 
0.0950 ± 0.0050 
0.0115 ±0.0005 
0.0306 ±0.0015 
0.0422 ±0.0058 
0.1341 ±0.0141 
0.2482 ±0.01 18 
0.1531 ±0.0081 
0.0502 ± 0.0047 
0.0423 ± 0.0043 
0.1224 ±0.0064 
0.1351 ±0.0071 
0.0241 ±0.0012 
0.0301 ±0.0015 
0.1910 ±0.0196 
0.0057 ± 0.0003 
0.0076 ± 0.0004 
0.0611 ±0.0029 
0.1339 ±0.0127 
0.0370 ±0.0018 
0.0909 ± 0.0091 
0.1139 ±0.0111 
0.0670 ± 0.0070 
0.1243 ±0.0062 
0.1339 ±0.0127 
0.2511 ±0.0239 
0.0456 ± 0.0024 
0.0832 ± 0.0040 
0.0882 ± 0.0042 
0.0905 ± 0.0043 
0.0906 ±0.0088 
0.0059 ± 0.0003 
0.0212 ±0.0022 
0.4022 ± 0.0422 
0.0766 ± 0.0075 
0.2009 ±0.0191 
0.0923 ± 0.0090 
0.0753 ± 0.0072 
0.1206±0.0115 



